On the interface polaron formation in organic field-effect transistors 
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A model describing the low density carrier state in an organic single crystal FET with high-K gate 
dielectrics is studied. The interplay between charge carrier coupling with inter-molecular vibrations 
in the bulk of the organic material and the long-range interaction induced at the interface with a 
polar dielectric is investigated. This interplay is responsible for the stabilization of a polaronic state 
with an internal structure extending on few lattice sites, at much lower coupling strengths than 
expected from the polar interaction alone. This effect could give rise to polaron self-trapping in 
high-K organic FET's without invoking unphysically large values of the carrier interface interaction. 
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I. INTRODUCTION 

Organic field-effect transistors (OFETs) are the ele- 
mentary building blocks of "plastic electronics" .^ In these 
devices, charges are induced in a thin conducting channel 
at the interface between an organic semiconductor and a 
gate insulator. While commercial applications have al- 
ready appeared, mostly based on lovif-cost organic thin- 
films, the fundamental physical mechanisms governing 
the charge dynamics in these devices are not fully un- 
derstood. As it has become clear in recent years, even 
when high purity crystalline semiconductors are used in 
order to minimize the effects of structural and chemi- 
cal disorder, the electronic characteristics of OFETs are 
strongly affected by interactions taking place both within 
the semiconducting material as well as with its close en- 
vironment. This can be traced back to the extremely 
narrow electronic bands resulting from the weak van der 
Waals intermolecular bonding, that make organic semi- 
conductors much more sensitive to electronic interactions 
than their inorganic counterparts. For this reason, pola- 
ronic effects have been shown to play an important role 
in these devices. 

Several electron-lattice interaction mechanisms rele- 
vant to organic molecular crystals have been identified 
and studied in the past. These include Holstein-like cou- 
plings to intra-molecular vibrations, leading to a local 
reorganization of the molecular energy levels^r— as well 
as Su-Schrieffer-Heeger (SSH) or Peierls-like couplings 
where the transfer integrals between the molecules are 
modulated by the inter-molecular motiouj^^— The lat- 
ter mechanism couples directly the electronic motion to 
the strong lattice fiuctuations arising from the mechan- 
ical softness of these compounds, and has been recently 
identified as the dominant limiting factor of the charge 
mobility in pure crystalline semiconductors i^SrJi 

In addition to such intrinsic mechanisms, interface- 
related effects have also been demonstrated in several 
systematic studies of organic FETs using different gate 
dielectrics. Sources of charge trapping have been identi- 
fied, either related to the interface qualit y ^^d^ or to the 
long-range forces arising from polar impurities at the or- 
ganic/dielectric interfaceJi"— Furthermore, it has been 



shown that the long-range polarization induced in high-K 
gate materials by the charge carriers themselves can lead 
to self-trapping of the carriers. Even when sufficient care 
is taken in order to minimize extrinsic sources of disorder, 
Frohlich polarons are formed due to the polar long-range 
interaction of the carriers in the organic semiconductor 
with the phonons at the semiconductor/dielectric inter- 
face. The strength of this remote electron-phonon (e- 
ph) interaction can be tuned by an appropriate choice 
of the gate dielectric. The "metallic-like" carrier mobil- 
ity characteristic of crystalline organic semiconductors 
can be converted into a much lower, thermally activated 
mobihtyi^"— 

While there is a vast theoretical literature dedicated 
to the different electron-lattice interaction mechanisms 
mentioned above, up to now there have been no attempts 
to study the interplay between bulk and interface e-ph 
couplings. To fill this gap, in this work we analyze a 
model that includes both bulk and interface effects as can 
be realized in organic FETs, and treat them in a common 
framework. We focus on the combined effects of the SSH 
coupling of the electronic transfer integrals to the inter- 
molecular vibrations and of a Frohlich long-range e-ph 
interaction as induced by the presence of the gate dielec- 
tric in a field-effect geometry. Apart from its relevance to 
organic/dielectric interfaces, from the fundamental view- 
point this model presents an interesting combination of 
two o priori competitive mechanisms. Our results show 
that a rather weak SSH coupling strength as estimated 
in Refs. lldllH can have an unexpectedly strong effect in 
stabilizing a polaronic state when a moderate long-range 
polar couphng is present. Therefore, self-trapped states 
of small radius can exist at common organic/dielectric 
interfaces, even in such cases where the carrier-interface 
phonon interaction alone would not be sufficient to pro- 
duce polaronic self-localizatior^. This provides a mi- 
croscopic theoretical basis for the experimental results 
of Ref. [2l|, where a finite activation energy indicative 
of self-trapped states was observed using gate materials 
such as AI2O3 and Ta205. 

This paper is organized as follows. The model under 
study is introduced in Section [Hi The two methods of 
solution that will be used are described in Section Hill In 



Section IIVI we introduce the main quantities of interest 
and present the resuhs of our calculations. Section FVl is 
devoted to the conclusions. A detailed derivation of the 
form of the long-range interaction with the gate material 
is presented in Appendix [^ 



II. THE MODEL 

We consider the following one dimensional tight- 
binding model 

H ^ -iX^t-*- + "("^ + "^ " "'■+! ^ "l+i)] (4ci+i + H.c.) 
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where electrons or holes described by the creation and 
destruction operators c\ , Ci move on a lattice labelled by 
site index i. These interact with molecular displacements 
Xi — al + Gi via a SSH interaction, describing the trans- 
fer integral modulation on the distance between near- 
est neighbors with strength a, as well as with optical 
modes Yq ^ bq + b_q at the polar interface via a coupling 
Mq. To keep the discussion as simple as possible and il- 
lustrate the main consequences of the interplay between 
the "off-diagonal" short range SSH coupling and the "di- 
agonal" long-range polar coupling, we shall restrict the 
electronic motion to a one-dimensional chain. Yet, our 
results can be considered as being qualitatively represen- 
tative of organic/dielectric interfaces in the more realistic 
two-dimensional case. 

As it is discussed in Appendix A we consider the one 
dimensional interface model introduced in Refs. [24l | 25| . 
The e-ph interaction matrix element reads in this case 
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with i?o a cut-off length of the order of the lattice spac- 
ing a and q a strictly one-dimensional vector. As it is 
clear from Eq. ([2]), the interaction between the elec- 
tron charge and the associated lattice polarization is long 
ranged, with a power law decay in real space. We define 
two dimensionless coupling constants for the two different 
el-ph couplings: 



Ala — 



^ 2ujLB.t 



and 



A 



SSH 



^SSH 



(3) 



(4) 



We shall study the system in the range of parameters 
typical of single crystal organic semiconductors, taking 



rubrene as a case study. We set ujlr = 0.25i as an aver- 
age value for the interface phonons in the case of common 
polarizable gate materials^i and ujssh — 0.05t for the 
relevant intermolecular vibrations^^'^'* {Ti = 1 through- 
out the paper) . The lattice cutoff parameter entering in 
Eq. © is chosen to be Rq = O.Sa.^^ 

In view of the small values of lulu and ujssh, we will 
study the system by using two approaches: the adiabatic 
approximation, valid for vanishing ojlr and ujssh, and a 
variational exact diagonalization (VED) method able to 
include non adiabatic contributions. 



III. METHODS 
A. Adiabatic approximation 

For values of ujlr and ujssh that are much smaller 
than t, the lattice polarization cannot follow the elec- 
tronic oscillations and the wave function of the system 
can be factorized^^ into a product of normalized varia- 
tional functions \ip) and |/) depending on the electron 
and phonon coordinates respectively: 
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This factorization becomes exact when we neglect com- 
pletely the ionic kinetic energy i.e. ujssh, ujlr — >> . In 
Eq.®, |0)gj is the electron vacuum state and (j){Rm) are 
parameters that satisfy the relation: 
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The expectation value of the Hamiltonian ^ on the 
state \(p) provides an effective Hamiltonian for the lattice 
fields: 

{^\H\^) = -tY,W {Rrn)(t){R,n + 1) + h.c] 
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and 
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V^ = -at U{R,n+l) - HRrn-l)) (f>* (Rm) + h.c] . (10) 



It can be proven that the ground state of this Hamil- 
tonian is a product of two coherent states: 
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where |0)^^ and |0)gg^ are the vacuum states of the two 
lattice fields. We stress that the parameters Vm are real 
[seeEq.lfTOll]. 

At this stage, the phonon state |/) can be used to 
obtain an effective Hamiltonian for the charge carrier: 

-^e// = -i^[l + 2a(t), - Wi+i)] (c|Q+i+cJ^iCi) 

i 

- Y.^UY.^M,e^''''Pl + h.c.)+K (12) 
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where 
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The present adiabatic approximation breaks transla- 
tional invariance. The problem is reduced to studying 
an electron, confined by the potential well raised by the 
LR interaction. The hopping amplitude from site to site 
within the well depends on the specific site pair, via the 
SSH coupling. The wavefunction parameters 4'{Rm) are 
determined variationally. 



B. The VED diagonalization technique 

In the systems on which we are focusing our attention 
the phonon frequencies lossh and ujlr are small com- 
pared to the transfer integralsi^i The adiabatic approach 
discussed above should therefore provide a reasonably ap- 
propriate description of the ground state properties. It is 
known, however, that a proper inclusion of non adiabatic 
terms in the calculation removes the fictitious transla- 
tional symmetry breaking, which is found within the pure 
adiabatic approximation, and can introduce significative 
changes at any non-zero phonon frequencies. Further- 
more physical quantities related to the quantum nature 
of the phonons (number of phonons in the ground state) 
or the delocalization of the electronic wave function (ef- 
fective mass and ground state spectral weight) can be 
accessed only via a non adiabatic approach. To this aim 
we implement here a method, combining a variational 
approach and an exact diagonalization technique, which 
is able to provide an accurate description of the ground 
state properties at all values of the system parameters. 
We shall first analyze the two interactions separately, and 
then introduce the method for treating the full Hamilto- 
nian Eq. (P. 

It has been shown22r^ that the ground state features 
in the case of "diagonal" electron-lattice interactions, 
such as the LR term Eq. ([2]), are very well reproduced 
by methods assuming as starting point the approach in- 
troduced by Toyozawa.^- The idea is to restore a trans- 
lationally invariant Bloch state using the wavefunctions 



provided by the adiabatic approximation, i.e. by taking 
a superposition of localized states centered on different 
lattice sites. The trial wave function that accounts for 
the translational symmetry reads: 
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In Eqs. (ll4liroi) the index n denotes the distance between 
the electron position and the lattice modes. 

The next step is to generalize this wavefunction to cor- 
rectly recover the weak coupling regime where the elec- 
tron is weakly dressed by lattice oscillations, as predicted 
by perturbation theory. This can be achieved by treating 
the parameters entering into the definition of the function 
Pq as independent variational parameters;^ rather than 
being constrained by the shape of the electron function as 
in Eq. ^. Actually, also in the perturbative regime the 
lattice displacement is described by a coherent state, but 
with a different function pq including an electron recoil 
term due to the scattering with phonons. ^^'^^ The pa- 
rameters (j)k{Ri) have to be determined via numerical di- 
agonalization of the corresponding Hamiltonian, defining 
the Variational Exact Diagonalization method described 
below. In this way we obtain a wavefunction that is able 
to capture the weak and strong coupling limits and in- 
terpolates between them in the intermediate regime. 

In practice we use a basis where each vector is a lin- 
ear superposition, with appropriate phase factors, of the 
translational copies of a state having a given electron po- 
sition and lattice configurations (phonons and electron 
are translated together, with their relative position be- 
ing held fixed). We choose as phonon states those ap- 
pearing in Eq. (|16p where pq is a variational function. 
There are then M basis vectors, labelled by the n in- 
dex appearing in Eq. (J16p . and each of them includes 
the variational function pq. The coefficients (j)k{Ri) in 
Eq. ([T5)) are determined for each fixed set of the parame- 
ters defining the variational function pq by diagonalizing 
the Hamiltonian matrix in the chosen basis. These do 
not coincide, in general, with the coefficients of the func- 
tion Pq in Eq.®. Summarizing, we have to diagonalize 
a MxM matrix for each value of the variational function 
Pq. Finally, in the following, we assume that the varia- 
tional parameters (j){Ri), entering into the definition of 
the function pq in Eq. Q, are different from zero up 
to 3rd nearest neighbors. We have checked that this as- 
sumption allows to obtain convergency for all the values 
investigated of the model parameters. 

The procedure described above could be straightfor- 
wardly generalized to include also the interaction with 



SSH phonons. Unfortunately, it can be shown that it 
does not capture the correct physical features of the SSH 
Hamiltonian in the absence of LR interaction.— To treat 
the SSH term we therefore apply a different approach 
based on an exact diagonalization procedure, where we 
do not assume as basis set the coherent states of the 
adiabatic approach. In this approach the phonon vari- 
ables are efficiently limited to a new set which is able 
to describe successfully the ground state properties, as it 
will be shown below. The real bottleneck conies from the 
phonon Hilbert space, that is infinite. To circumvent this 
difficulty, we keep states where the lattice deformations 
involve only the sites up to 3rd nearest neighbors of the 
lattice site where electron is located, all the others being 
undeformed 
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where n indicates the electron position. We apply a trun- 
cation procedure considering only the states with at most 
L quanta, i.e. J2j''^j — ^' where L is fixed by requir- 
ing convergency of the ground state energy. These states 
provide a correct physical description of the ground state 
properties, since the lattice deformations are effectively 
localized around the site where the electron is located. 
Furthermore, in analogy with the LR coupling, we use the 
translational symmetry associated to periodic boundary 
conditions to reduce the size of the sparse Hamiltonian 
matrix, requiring that the states have a definite momen- 
tum. 

When both LR and SSH couplings are considered, we 
adopt a VED diagonalization technique on a truncated 
phonon basis that combines the two subspaces introduced 
separately for the a = and g = cases. This ba- 
sis is given by the product of the vectors describing the 
two different couphngs separately, Eqs. p^ and (IT71) . 
The ground state is found by using the modified Lanc- 
zos algorithm. ■^^ We emphasize that this approach retains 
the full translational invariance requiring that each basis 
vector has a well defined momentum. 



IV. RESULTS 

Due to their different microscopic nature, the two in- 
teraction terms in Eq. (1) induce different types of po- 
laron states: the LR coupling to the electron density fa- 
vors localized states centered on the lattice sites (site po- 
laron), while the SSH coupling to the electronic transfer 
integrals tends to localize the polarons on inter-molecular 
bonds, i.e. with the electronic charge centered between 
two neighboring molecular sites (bond polaron). This 
competition results in the phase diagrams reported in 
Fig. [l] which constitute the central results of this paper. 
The most notable feature is the existence of a reentrant 
phase where a bond polaron state is stabilized by the 
interplay between the two microscopic interaction mech- 
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FIG. 1: Phase diagram a) in the adiabatic approximation and 
b) in the VED approximation. The dashed (blue) line marks 
the crossover from large (LP) to small (SP) site-centered po- 
laron probed by the bimodality of the LR-PDF. The solid 
(red) line indicates the crossover to the bond polaron (BP) 
state probed by the bimodality of the SSH-PDF while the 
dashed (green) line is the transition to the bond polaron state 
probed by the electronic wavefunction (only in panel a). In 
panel a) the dashed double dotted red line is the analytic 
estimate of Eq. ((20)) . 



anisms. This result is robust against the presence of 
phonon quantum fluctuations, as is apparent from the 
comparison of the phase diagram obtained from the adi- 
abatic approximation, valid in the limit of vanishing fre- 
quencies uJssH,^LR ^ [FiglT](a)], with the results at 
finite phonon frequencies [FiglT](b)]. 



A. Adiabatic approximation 

In the adiabatic approximation the ground state elec- 
tronic wave function is always localized. However, for 
small values of Xssh, a crossover from a large to a small 
site-centered polaron occurs upon increasing the LR in- 
teraction strength Xlr [blue dashed line in Fig. [T] (a)]: 
the wavefunction extension is progressively reduced, be- 
ing always centered on a given molecular site. The locus 
of the crossover can be determined by analyzing the qual- 
itative changes induced on the phonon wavefunction, as 
described at the end of this Section. The large to small 
polaron crossover is only weakly affected by the presence 
of the SSH interaction. 

As XssH increases, the symmetry of the electronic 
wave function changes continuously from site-centered to 
bond-centered. The breaking of translational invariance 
allows us to analyze this transition by direct inspection 
of the electronic wavefunction, which is shown in Fig. 
[2] together with its associated distortions in the bond- 
centered (left panel) and site-centered polaron (right 
panel) states. The transition between the two phases 
is indicated by the green dashed line in FiglTJi. 

Contrary to the large/small polaron crossover, that is 
rather insensitive to the value of Xssh , the site to bond- 
polaron transition is significantly modified by the pres- 
ence of the LR interaction, which strongly stabilizes the 



bond polaron state. The boundary of the bond-centered 
polaron phase shows a markedly reentrant behavior, 
reaching a minimum Xssh = 0.072 at Xlr = 0.4.^^ 
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FIG. 2: Bond polaron properties for Xssh = 0.4 Xlr = f .0 
(left panels) and Xlr = 1.5 (right panels) in the adiabatic 
approximation. Upper panel: the adiabatic electronic proba- 
bility. Central panel: SSH average distortions of the organic 
molecular lattice. Lower panel: LR average distortions at the 
organic/dielectric interface. 



It is instructive to understand the stabilization of the 
bond polaron phase induced by the long-range interface 
interactions from the following energy balance argument. 
In the absence of intermolecular coupling {Xssh = 0) the 
energy of the large polaron can be written in general as 



Elp = -2t - AE 



LP 



(18) 



where the second term can be expanded for small Xlr 
as AElp = t{aiXLR + a2A|^ -I- ...). On the other hand 
we see from Fig. [5] that in the bond polaron phase the 
adiabatic electronic wave function is mainly localized on 
two neighboring sites, with exponentially vanishing tails 
on the adjacent molecules. Neglecting such tails provides 
the following estimate for the bond polaron energy: 
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where the two contributions come respectively from 
SSH and LR interactions. The parameter rj == 



EjM,Pe»VE, 



\M,\' 



{a the lattice spacing, < 77 < 



1) takes into account the spatial variation of the long- 
range potential- well between two neighboring sites. Com- 
paring the two energies we obtain a transition when 
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and a reentrant behavior is found when ai < (1 + r/)/2. 
Observing that rj is largest for slowly varying potentials, 
we recognize that the reentrance is strongly favored by 
the long range nature of the polar dielectric interaction. 
The reason is that such long-ranged polar interaction, 
for which the preferred polaronic state would be site- 
centered, is not extremely sensitive to the local details 
of the wavefunction, so that it can still provide a sizable 
energy stabilization in the case of a dimer structure. It 



should be stressed that the arguments presented here are 
only qualitative. For a quantitative analysis of the large 
polaron to bond polaron phase boundary, the extended 
nature of the bond polaron wavefunction must be taken 
into account, which further enhances the reentrant be- 
havior as compared to Eq. ([20| (see FigH]). 

Phonon related characteristics also change across the 
bond polaron transition, following the changes of the elec- 
tronic wavefunction. As can be seen from Fig. [51 the 
average SSH distortions within a bond polaron are op- 
posite, since this configuration minimizes the SSH inter- 
action energy. Within the entire bond polaron region of 
the phase diagram, even at the relatively large value of 
the LR couphng Xlr = 1 illustrated in Fig. [2l the LR 
distortions follow the charge distribution determined by 
minimizing the SSH interaction energy. 

The phonon probability distribution functions (PDF) 
provide an alternative tool to characterize the ground 
state propertieSji^ with the advantage of being directly 
generalizable to the translationally invariant ground 
states examined in the next section. The SSH-PDF, 
which represents the distribution of SSH displacements 
at site i + j given an electron on site i , is defined as 
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where |Xi-(-j > is the coordinate eigenstate correspond- 
ing to the displacement of SSH phonons on site i + j. 
The LR-PDF is defined by substituting X with Y in Eq. 
(PT|) . It is worth noting that the PDF is closely related to 
the distributions of interatomic distances that can be di- 
rectly measured by neutron spectroscopy such as pulsed 
neutron diffractioui^ 

The local {j = 0) part of these functions is plotted 
in Fig. [3] for different polaronic ground states. In the 
small Xssh region of the phase diagram [see Fig. [Da)], 
the large-to-small polaron crossover is highlighted by the 
bimodality in the LR-PDF [panel (a) of Fig. |3]. On the 
other hand, when the bond polaron state is well formed, 
the SSH-PDF shows two symmetric peaks [panel (b)], 
whose distance increases with Xssh- This provides an 
alternative method to locate the bond polaron region, 
which is shown in Fig. [l](a) as a red continuous line. The 
resulting phase diagram shows that a broad crossover 
region emerges (between the green and red lines in Fig. [T] 
where the electronic wavefunction has already undergone 
the bond transition but a real bimodality of the SSH- 
PDF is still not observed, with the SSH-PDF showing 
large non-Gaussian fiuctuations. Let us note that in the 
bond polaron region of the phase diagram, the analysis 
of the LR-PDF does not show any evidence of a large 
to small bond polaron crossover, at least for values of 
Xssh < 0.5. 

Summarizing, the adiabatic phase diagram shows three 
distinct regimes: small site-centered polaron, large site- 
centered polaron and bond polaron. The bond polaron 
regime, with the electron wavefunction mostly localized 





FIG. 3: (left) Local probability distribution function Pji^JiiY) 
of the LR phonons (see text) at Xssh = 0, Xlr = 0.7 (blue 
dotted line) and Xlr ~ 1.5 (red solid line); (right) Same for 

the SSH phonons: Ps^sni^) at Xssh = 0.4, Xlr = (blue 
dotted line) and Xlr = 1 (red solid line) 



In Fig. Hwe plot the SSH-PDF and LR-PDF for j = 

— 1,0, 1, in the small Xssh region (we use Xssh = 0.2) 
of the phase diagram, across the large-to-small polaron 
crossover [see Fig. [I](b)]. This crossover is evinced by 
the bimodality in the LR-PDF, occurring at Xlr — 1.4. 
Both the on-site {j = 0) and the inter-site LR-PDF {j — 

— 1,1, which are equal by symmetry), become bimodal at 
the crossover. In this region, the SSH phonons present 
a nearly Gaussian distribution with opposite distortions 
around the polaron localization site. It is interesting to 
note that the average value of the SSH displacement has 
a non-monotonic dependence on the LR interaction and 
it is largest for intermediate values of Xlb., signaling the 
proximity to the bond polaron regime. 



on a molecular dimer, is stabilized by the presence of 
even a moderate long range coupling. 



B. Variational Exact Diagonalization method 

Fig. [T] (b) illustrates the phase diagram obtained us- 
ing the VED scheme described in section IIIIBI Trans- 
lational symmetry breaking is prevented by the quan- 
tum fluctuations of both SSH and LR phonons. Accord- 
ingly, all the lines shown in Fig. [1] (b) correspond to 
crossover lines separating translationally invariant states 
with different characteristics, as obtained by analyzing 
the phonon PDFs.^ 
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FIG. 4: SSH (left) and LR (right): phonon PDF at Xssh = 
0.2: upper panel: Xlr = 1.0, central panel: Xlr ~ 1.5, 
lower panel: Xlr = 2.5. Solid (red) line represents j = (see 
text); dashed blue (double dashed black) line represents j — I 

U = -1). 
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FIG. 5: SSH (left) and LR (right) phonon PDF at Xssh = 0.4 
in the bond polaron regime. The solid (red) line represents 
the on-site contribution, j = (see text); dashed blue (double 
dashed black) line represents the inter-site part j — 1 {j = 

In the bond polaron phase both the LR and the SSH- 
PDF have a bimodal character, as shown in Fig. \5\ The 
SSH-PDF displays bimodality in both the on-site (Pssh) 
and nearest neighbor (Pssh) contributions. The site 
where the polaron is located belongs to either a left or 
right-centered dimer, a situation which is consistent with 
the translational invariance of the ground state. The 
LR-PDF instead is bimodal only on the nearest neighbor 
{P^R )j with two peaks corresponding to distorted and 
non-distorted sites respectively. The presence of a distor- 
tion with the same amplitude on nearest neighboring sites 
is another signal of the tendency to form a dimer. This 
tendency is also displayed by the LR distortions which, 
in this region of the phase space, follow the charge dis- 
tribution dictated by the SSH interaction. 

Let us note that the reentrant bond polaron phase at 
^LR 7^ has a different symmetry than that at Xlr = 0. 
Indeed, according to Ref.|37|, the transition to the bond 
polaron phase induces the ground state wavenumber to 
change continuously from k = to k = ±7r/2 in a nar- 
row region around Xssh = 0.5. As soon as Xlr is dif- 
ferent from zero, we do not find any evidence of a fc 7^ 
ground state for Xssh < 0.5. To illustrate this point, 
the electron dispersion characteristic of this phase, where 
the minimum is located at fc = 0, is plotted in Fig. [SI 
We conclude that in the reentrant bond polaron phase, 
the LR coupling stabilizes a Bloch wavefunction having 
k = 0. 

In Fig. [7] we plot the inverse effective mass and the 





FIG. 6: Polaron band at Ala = 1.0 and Xssh 
(reentrant bond polaron phase), in units off. 
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0.288 FIG. 8: Electron kinetic energy, in units of the bare kinetic 
energy {Xssh ~ Xlr = 0), as a function of \lr for different 
values of Xssh- 



ergy stays relatively constant over a broad range of Xlr- 
This plateau is another clear signature of the bond po- 
laron phase. Here the effective mass and the quasiparticle 
spectral weight are very small, indicative of polaron self- 
trapping, yet a large contribution to the kinetic energy 
arises from the internal dimer structure of the polaron, 
with the electron moving freely along the intermolecular 
localization bond. 



FIG. 7: Inverse effective mass renormalization (upper panel, 
m is the band mass) and the quasiparticle spectral weight 
(lower panel) as a function of Xlr for different values of Xssh- 



quasiparticle weight, defined as Z = {GS\ cJ,^q|0 
a function of Xlb. for different values of the SSH coupling 
^SSH- At small vaues of Xssh the inverse effective mass 
first decreases smoothly with Xlr and then drops at the 
large-to-small polaron crossover around Xlr ~ 1.4— 1.5. 
The quasiparticle spectral weight shows a similar depen- 
dence on Xlr, although with a steeper reduction which is 
characteristic of non-local e-ph interactions j^^*^ In con- 
trast with the rather smooth Xlr dependence, a much 
stronger dependence on Xssh can be seen in Fig. [71 
signaling the crossover into the reentrant bond polaron 
regime at Xssh < 0.5. We note that the quasiparticle 
spectral weight and inverse effective mass have in this 
case the same behavior with Xlr as it occurs for local 
e-ph interactions;^ i.e. they tend to be mutually pro- 
portional. 

Complementary insight on the electron dynam- 
ics in the different phases can be obtained from 
the analysis of the total kinetic energy, Ekin = 

-tY,i \{4^i+i + 4+i^i))' shown in Fig. 13 Unlike the 
quasiparticle spectral weight Z and effective mass to*, 
this quantity gets contributions from both the coherent 
long distance motion of the polaron, and from the inter- 
nal motion of the electron within the lattice potential- 
well. The large-to-small polaron crossover at small Xssh 
causes a sharp drop of Ekin around Xlr — 1.5. For 
Xssh > 0.25 a region emerges where the kinetic en- 
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FIG. 9: Average number of SSH (upper panel) and LR 
phonons (lower panel) in the polaronic ground state as func- 
tion of Xlr for several values of Xssh- 

The large number of SSH phonons participating in the 
polaronic state in the reentrant phase, shown in Fig. 
[21 indicates that a strong coupling with intermolecular 
modes in the organic crystal can be effectively triggered 
even at relatively weak values of Xssh, by suitably turn- 
ing on the LR interaction. At the same time, the number 
of LR phonons stays relatively small, exhibiting only a 
weak dependence on Xssh (see lower panel of Fig. [S]). 



V. CONCLUSIONS 

In this paper, we have investigated a microscopic 
model that describes the interplay between the inter- 
molecular electron-phonon coupling and the additional 
interaction arising in crystalline organic semiconductors 
due to the long range polarization at the interface with a 
gate dielectric, two ingredients that are commonly found 
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in organic field effect transistors. The ground state prop- 
erties obtained by two complementary approaches — 
the adiabatic approximation and a more refined varia- 
tional method that correctly preserves translational in- 
variance in the presence of phonon quantum fluctuations 
— demonstrate the synergistic cooperation between bulk 
and interfacial interactions. In particular, we have found 
that a large-mass polaronic state with a dimer structure 
and strong inter-molecular distance fluctuations is stabi- 
lized at relatively low values of both the long-range po- 
lar coupling Aii{ and the intermolecular coupling Xssh- 
Such bond polaron reentrant phase arises in a region of 
the phase diagram where the two interactions taken sep- 
arately would not give rise to strong polaronic effects. 
The long range nature of the polar interaction is cru- 
cial for the establishment of this cooperative interplay. 
The latter occurs because the long range interface polar- 
ization is able to provide a sizable energy gain without 
much perturbing the local dimer structure dictated by 
the SSH intermolecular coupling, so that the system can 
take advantage of both microscopic interaction mecha- 
nisms simultaneously. 

From the device perspective, our findings provide a 
microscopic basis to explain the existence of small-radius 
self-trapped states without invoking unrealistic values of 
the polar coupling. Indeed, the electrical characteristics 
of high-K organic FETs have been analyzed in Ref.j21l] 
in terms of the polar interface interaction alone, and the 
observed thermally activated mobilities are not entirely 
consistent with the estimated values of Xlr, as has been 
pointed out in Ref.[19|. The present results show that 
the proper inclusion of a bulk electron-lattice interaction 
mechanism such as the intermolecular SSH coupling ex- 
tends the region of small-radius polaronic self-trapping 
down to much lower values of \lr- Let us stress that 
an analogous interplay as the one studied here can be 
expected when combining the effects of a bulk inter- 
molecular electron-lattice coupling and the interaction 
with localized polarizable impurities such as those con- 
sidered in Refs. [3il. 

Comparing the inter-molecular coupling strength esti- 
mated for Rubrene, Xssh — 0.17fii with the phase di- 
agrams of Fig. [1] we conclude that interfaces involving 
this organic material and the aforementioned gate mate- 
rials lie in a region of the phase diagram where the in- 
termolecular/gate synergistic effects are expected to be 
of relevance. More importantly, our results could ex- 
plain why at present there are only few organic semi- 
conductors (among which, Rubrene) where metallic-like 
mobilities have been measured in a FET configuration. 
Other materials with slightly larger values of the inter- 
molecular coupling strength, Xssh, would rapidly enter 
the bond-polaron phase, causing the carriers to undergo a 
self-trapping transition as soon as a moderate interaction 
with the dielectric environment is present. 

Although the model investigated in this paper includes 
some drastic approximations, i.e. the model is one- 
dimensional and the real lattice structure of the organic 



FET is not taken into account, we believe that the main 
conclusions provided here will not be qualitatively mod- 
ified by the inclusion of more realistic lattice structures 
and interactions. The detailed study of the transport 
characteristics in the bond polaron phase identified here 
remains as an open question for future work. 



Appendix A: Long-range interaction at the interface. 

In this work we studied the interplay between SSH 
and Frohlich e-ph interactions in one-dimension. While 
the definition of the SSH coupling in one-dimension is 
unambiguous,^^ the polar e-ph interaction term depends 
on the geometry of the interface. We now examine two 
physically different cases. 

The case that is studied in the paper is a truly 
one-dimensional interface, derived for a wire of square 
cross-section and permittivity ei embedded in a three- 
dimensional polar dielectric with permittivity €2 .—1^ We 
specialize the results of Ref.[25] by allowing the motion 
of the charges only on the symmetry axis of the interface 
z with X — 0,y = 0, atsi distance Rq from the polar 
dielectric with permittivity 62- The medium with per- 
mittivity ei is taken to be the vacuum. In this geometry 
the interaction of electrons with the interface phonons 
reads 



M„ = Mqc 



-qRo 



with Rq = {\x\ + \y\)/^/2 or, in real space, 

Rq 



/(i?)(X 



i?2 + i?2 



(Al) 



(A2) 



Notice that here the phononic modes are purely one 
dimensional modes, i.e. they are propagating waves only 
along the wire (direction z). 

It is interesting to compare the one dimensional 
interaction discussed above with a two dimensional 
model. The model that we consider is the generaliza- 
tion of the case studied previously by Mori and Andoj^ 
that describes the coupling between electrons and two- 
dimensional interface phonons. In this case the inter- 
action matrix element in momentum space reads Mq — 
Moe~''^ /y/q, with the distance z between the electrons 
and the interface acting as a short-distance cutoff. Per- 
forming the Fourier transform we find the Hypergeomet- 
ric function 



fiR) 



MqTT 



3/2 



z3/2 



-2Fi(3/4,5/4,l,-i?7z2) 



(A3) 



This scales as 1/ R^/^ for large R ^ z, and tends to a 
constant /(O) = Mq^/ z'^/'^ for R<^ z. This function is 
well approximated by the power-law form 



I{R) = 



M'^ 



(i?2 + i?2)3/4 



(A4) 



where Rq = {Ay/Trc)'^/^z ~ 0.66z and M(^ = cMq. It can 
be easily checked by retransforming back to k-space that 
this power-law form gives rise to the correct q —> limit 
and lies within 10% of the original one for all q ^ I/Rq. 
Then, in both one and two-dimensional cases, Eqs. 
()A2|A4p give rise to long range interactions, describing 
qualitatively similar physical situations. We have chosen 
to study explicitly the first model in the body of the pa- 



per, as in that case the properties of the phononic degrees 
of freedom can be directly accessed by the solution of a 
genuinely one dimensional model. 
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